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Abstract 

We summarize the current view on Parton-Hadron duality as it applies to B meson 
decays. It is emphasized that an OPE treatment is essential for properly formulating 
duality and its limitations. Duality violations are unlikely to become the limiting 
factor in describing semileptonic B width vis-a-vie higher order corrections. The 
consistent extraction of the b quark mass from B production and decays provides a 
striking example of the theoretical control achieved. 



^On leave from Physics Dept., University of Notre Dame du Lac, Notre Dame, IN 46556, U.S.A. 



1 Introduction 



Parton-hadron duality | - or duality for short - is one of the central concepts in contem- 
porary particle physics. It is invoked to connect quantities evaluated on the quark-gluon 
level to the (observable) world of hadrons. It is used all the time, more often than not 
without explicit reference to it. A striking example of the confidence the HEP commu- 
nity has in the asymptotic validity of duality was provided by the discussion of the width 
r(Z° — > HbH'ijX). There was about a 2% difference in the predicted and measured decay 
width, which lead to lively debates on its significance vis-a-vis the experimental error. No 
concern was expressed about the fact that the width was calculated on the quark- 
gluon level, yet measured for hadrons. Likewise the strong coupling as{Mz) is routinely 
extracted from the perturbatively computed hadronic width with a stated theoretical 
uncertainty of 0.003 with translates into a theoretical error in ThadiZ'^) of about 0.1%. 

There are, however, several different versions and implementations of the concept of 
duality. The problem with invoking duality implicitly is that it is very often unclear which 
version is used. In B physics - in particular when determining |V^(c6)| and - the 

measurements have become so precise that theory can no longer hide behind experimental 
errors. To estimate theoretical uncertainties in a meaningful way one has to give clear 
meaning to the concept of duality; only then can one analyze its limitations. 

In response to the demands of B physics a considerable literature has been created 
on duality over the last few years, which we want to summarize. We will emphasize and 
illustrate the underlying principles; technical details can be found in the references we list. 

Duality for processes involving time-like momenta was first addressed theoretically in 
the late '70's in references and 0. We sketch here the argument of Poggio, Quinn and 
Weinberg since it contains several of the relevant elements in a nutshell. The cross section 
for e^e^ hadrons can be expressed through an operator product expansion (OPE) of 
two hadronic currents. One might be tempted to think that by invoking QCD's asymptotic 
freedom one can compute a{e~^e~ — ^ hadrons) for large cm. energies s ^ ^qcd in terms 
of quarks (and gluons) since it is shaped by short distance dynamics. However production 
thresholds like for charm induce singularities that vitiate such a straightforward computa- 
tion. This complication can be handled in the following way. Consider the correlator of 
two electromagnetic currents: 



where n(g^) can be written in terms of a spectral function p(s) using an unsubtracted 
dispersion relation: 



It is well known that, to leading order in aem, p{s) is related to the total cross section for 

^This name might be more appropriate than the more frequently used gwarfc-hadron duality since gluonic 
effects have to be included as well into the theoretical expressions. 



rf^xe^^^ {0\T{J,{x)MO))\Q) = M.-M')n(g') 



(1) 




(2) 
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e~^e hadrons 

ais) = ^pis) (3) 

Relying on QCD's asymptotic freedom one computes the correlator (|l]) in terms of quarks 
and gluons for s in the deep Euclidean domain |s| ^ ^qcd] s is chosen Euclidean so 
that one avoids a proximity to singularities induced by hadronic thresholds like for charm 
production etc. From the spectral function p calculated in the Euclidean regime one can 
infer the cross section for physical, namely Minkowskian values of s. However, one cannot 
obtain it as a point-for-point function of s, only averaged - or 'smeared' - over an energy 
interval, which can be written symbolically as follows: 



pisEuclid) / dsa{e^e~ hadrons) (4) 

J so 



This feature is immediately obvious: for the smooth s dependence p has to be compared 
to the measured cross section e+e" — >■ hadrons as a function of s, which has pronounced 
structures, in particular close to thresholds for cc- and 66-production. 

This simple illustration already points to the salient elements and features of duality 
and its limitations: 

• An OPE description in terms of quark and gluon degrees of freedom for the observable 
under study is required. 

• This OPE has to be constructed in the Euclidean domain. 

• Its results are analytically continued to the Minkowskian domain with the help of a 
dispersion relation. 

• This extrapolation implies some loss of information; i.e. in the notation given above 

/rphadronic\ ^ /rppartonic\ f r\ 

where {...)w denotes the smearing which is an average using a smooth weigth function 
w{s)\ it generalizes the simplistic use of a fixed energy interval: 

(...)^ = / ds...w{s) (6) 



• Some contributions that are quite insignificant in the Euclidean regime and there- 
fore cannot be captured through the OPE can become relevant after the analytical 
continuation to the Minkowskian domain, as explained later on. For that reason we 
have used the approximate rather than the equality sign in Eq.(^). 

• One can make few universal statements on the numerical validity of duality. How 
much and what kind of smearing is required depends on the specifics of the reaction 
under study. 
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The last item needs expanding right away. The degree to which (T^^''*""*'^)^ can be trusted 
as a theoretical description of the observable (T^°'^''°"'*'^)«, depends on the weight function, 
in particular its width. It can be broad compared to the structures that may appear 
in the hadronic spectral function, or it could be quite narrow, as an extreme case even 
w{s) ~ 6{s — So). It has become popular to refer to the first and second scenarios as global 
and local duality, respectively. Other authors use different names, and one can argue that 
this nomenclature is actually misleading. Below we will describe these items in more detail 
without attempting to impose ex cathedra a uniform nomenclature ourselves. 

Irrespective of names, a fundamental distinction concerning duality is often drawn be- 
tween semileptonic and nonleptonic widths. Since the former neccessarily involves smearing 
with a smooth weight function due to the integration over neutrino momenta, it is often 
argued that predictions for the former are fundamentally more trustworthy than for the 
latter. However, as we shall see, such a categorical distinction is overstated and artificial; 
also it is not needed for the discussion in the following chapters. Of much more relevance 
is the differentiation between distributions and fully integrated rates. 

No real progress beyond the more qualitative arguments of Refs. [|I| and [0 occurred 
for many years. For as long as one has very limited control over nonperturbative effects, 
there is little meaningful that can be said about duality violations. Yet this has changed 
for heavy flavour physics with the development of heavy quark expansions, since within 
this OPE framework we can assess nonperturbative effects as well as duality violations. 

The remainder of this note will be organized as follows: In the next section we shall 
give a more precise definition of what is meant with "Parton Hadron Duality" , which then 
allows us a discussion of possible violations of duality in section ^ Based on this we give 
hints on how to check the concept of duality in section ^, before presenting conclusions. 



2 What is Parton— Hadron Duality? 

In order to discuss possible violations of duality we have to give first a more precise defi- 
nition of this notion, which requires the introduction of some theoretical tools. We follow 
closely the arguments given in the extensive reviews of Ref. and The central 

ingredient into this definition is the method of the Wilsonian Operator Product Expansion 
(OPE) frequently used in field theory to perform a separation of scales. In practical terms 
this means that we can write 

/ d'xe^^^ {A\T {nx)rm \A) ^ Y: (^) CW'; X){A\OM)x (7) 

for = — oo. The following notation has been used: \A) denotes a state that could 
be the vacuum - as for e~^e~ —>■ hadrons considered above - or a i? meson when describing 
semileptonic beauty decays. denote electromagnetic and weak current operators for the 

^It can be noted that even the authors of Ref. and [Q - although very close in the substance as well 
as the spirit of their discussion - do not use exactly the same terminology concerning different aspects of 
duality. 
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former and the latter processes, respectively; for other decays like nonleptonic or radiative 
ones one employs different AS = 1 operators; the On are local operators of increasing 
dimension. The operator of lowest dimension yields the leading contribution. In e~^e~ 
annihilation it is the unit operator Oq = 1, for B decays Oq = bb. They produce (among 
other things) the naive partonic results. Yet the OPE allows us to systematically improve 
the naive partonic result. The coefficients c^*" contain the contributions from short distance 
dynamics calculated perturbatively based on QCD's asymptotic freedom. Following Wil- 
son's prescription a mass scale A has been introduced to separate long and short distance 
dynamics; both the coefficients and the matrix elements depend on it, their product of 
course not. 

The perturbative expansion takes the form 



and is performed in terms of quarks and gluons. The expectation values for the local 
operators provide the gateways through which nonperturbative dynamics enters. 

The crucial point is that the OPE result is obtained in the Euclidean domain and has to 
be continued analytically into the Minkowskian regime relating the OPE result to observ- 
able hadronic quantities. As long as QCD is the theory of the strong interactions, it does 
not exhibit unphysical singularities in the complex plane, and the analytical continu- 
ation will not induce additional contributions. To conclude: duality between (T^^'^''""*'^)ty 
and (T^"''*°"*'^)„, arises due to the existence of an OPE that is continued analytically. It 
is just a restatement of QCD's basic tenet as the theory of the strong interactions that 
hadronic observables can be expressed in terms of quark-gluon degrees of freedom provided 
all possible sources of corrections to the simple parton picture are properly accounted for. 
It is thus misleading to refer to duality as an additional assumption. 

Up to this point our discussion was quite generic. To specify it for semileptonic B 
decays one chooses the current to be the weak charged current driven by 6 — > c or 
b ^ u. The expansion parameter for inclusive semileptonic decays is given by the energy 
release ~ I /{mi, — rUc) [l/mh\ for 6 — c [6 — > w]. For the exclusive mode B — > IvD* it is 
l/mfe and l/rric with the latter yielding the numerically leading contributions. 



3 Duality Violations and Analytic Continuation 



One of the main applications of the heavy quark expansion is the reliable extraction of 
|y(c6)| and One wants to be able to arrive at a meaningful estimate of the theo- 

retical uncertainty in the values obtained. There are three obvious sources of theoretical 
errors: 

1. unknown terms of higher order in as; 

2. unknown terms of higher order in l/mq] 




(8) 
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3. uncertainties in the input parameters as, mq and the expectation values. 

Duahty violations constitute uncertainties over and above these; i.e. they represent contri- 
butions not accounted for due to 

• truncating these expansions at finite order and 

• limitations in the algorithm employed. 

These two effects are not unrelated. The first one means that the OPE in practice is 
insensitive to contributions of the type e~"*'^/ ^ with n denoting some hadronic scale; the 
second one reflects the fact that under an analytic continuation the term e"™'^/'^, which is 
quite irrelevant for Q = b - though not neccessarily for Q — c I - turns into an oscillating 
rather than suppressed term sin(mQ//x). 

Of course we do not have (yet) a full theory for duality and its violations. Yet we know 
that without an OPE the question of duahty is ill-posed. Furthermore in the last few years 
we have moved beyond the stage, where we could merely point to folklore. This progress 
has come about for the following reasons: 

• We have refined our understanding of the physical origins of duality violations as due 
to 

— hadronic thresholds; 

— so-called 'distant cuts'; 

— the suspect validity of l/rric expansions. 

• We understand the mathematical portals through which duality violations can enter, 
namely that the innocuous Euclidean quantity e~'^Ql^ transmogrifies itself into the 
much more virulent Minkowskian quantity sin(m(3//i) under analytical continuation. 

• The quantity e^™'?/^ is actually innocuous for beauty, yet not necessarily for charm 
quarks. The 'Euclidean' quantity Fd*(0) - the formfactor for B — > IvD"" at zero 
recoil -, which is given by an expansion in 1/mc, could be vulnerable to such a 
duality violation. However, the heavy mass expansion for exclusive quantities such 
as Fd*(0) is not directly given by an OPE, thus this argument may not apply in this 
case. 

• We have come up with an increasing array of field-theoretical toy models, chief among 
them the 't Hooft model, which is QCD in 1-1-1 dimensions in the limit of ~^ oo. 
It is solvable and thus allows an unequivocal comparison of the OPE result with the 
exact solution. 

• For the analysis of 6 ^ c transitions we also have the small- velocity expansion as a 
powerful tool. 
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We will not go into any details here, since they can be found in the literature The 
models do exhibit duality violations, but only highly suppressed ones conforming to general 
expectations. There had been claims of sizeable duality violations in the previous literature; 
those have been analyzed carefully in Ref. [Q, where their flaws are pointed out explicitly. 

Based on general expectations as well as on analyzing the models one finds that indeed 
duality violations are described by highly power suppressed 'oscillating' terms of the form 

T{mQ) ~ (^-^^ sin(mQA) (9) 

for some integer power k. More generally one can state: 

• The primary criterion for addressing duality violation is the existence of an OPE for 
the particular observable. 

• Duality will not be exact at finite masses. It represents an approximation the ac- 
curacy of which will increase with the energy scales in a way that depends on the 
process in question. 

• Limitations to duality can enter only in the form of an oscillating function of energy 
or mg (or have to be exponentially suppressed). Duality violations cannot be blamed 
for a systematic excess or deficit in the decay rates. For example, no duality violation 
can convert niq into Mhq in the full width parametrically, only for discrete values of 

• The OPE equally applies to semileptonic as well as nonleptonic decay rates. Like- 
wise both widths are subject to duality violations. The difference here is quantita- 
tive rather than qualitative; at finite heavy quark masses corrections are generally 
expected to be larger in the nonleptonic widths. In particular, duality violations 
there can be boosted by the accidental nearby presence of a narrow hadronic reso- 
nance. Similar effects could arise in semileptonic rates, but are expected to be highly 
suppressed there. 

• It is not necessary to have a proliferation of decay channels to reach the onset of 
duality, either approximate or asymptotic. Instructive examples are provided by the 
so-called small- velocity kinematics in semileptonic decays and by nonleptonic rates 
in the 't Hooft model. 

Putting everything together it has been estimated with considerable confidence - at 
least by the authors of Ref . |^ - that duality violations in the integrated semileptonic width 
of B mesons cannot exceed the fraction of a percent level. As such we do not envision it 
to ever become the limiting factor in extracting |l^(c6)| and since the uncertainties 

in the expression for the semileptonic width due to fixed higher order contributions will 
remain larger than this level. The oscillatory nature of duality violating contributions is 
a main ingredient in this conclusion. It also shows that duality violations could become 
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quite sizeable if an only partially integrated width - let alone a distribution - is considered. 
Generally, for distributions the expansion parameter is not the heavy mass, rather it is a 
quantity such as l/[mQ(l — x)] where x is e.g. the rescaled charged lepton energy of 
a semileptonic decay. From equation @ one would expect that contributions the form 
sin(mQ[l — x])/[mQ{l — x)]'' appear in differential distributions. 

4 How can we check the vahdity of Parton— Hadron 
Duality? 

If in the future we were to find a discrepancy between the measured and predicted values 
for, say, a CP asymmetry in B decays, we had to check very diligently all ingredients upon 
which the prediction was based, in particular the values for V{cb) and V{ub), before we 
could make a credible claim to have uncovered New Physics. This means one needs a 
measure for potential duality violations that is not based purely on theoretical arguments. 

Most theoretical uncertainties are systematical rather than statistical. As it is the case 
for experimental systematics the most convincing way to establish control over them is to 
determine the same quantity in independent ways and analyze their consistency. The heavy 
quark expansions lend themselves naturally to such an approach since it is their hallmark 
that they allow the description of numerous decay rates in terms of a handful of basic 
parameters, namely quark masses and hadronic expectation values. Again the situation 
is very similar as for the perturbative series: once the coupling constant is determined 
(e.g. as{Mz)) from a measurement, one may use this as an input to all other perturbative 
calculations, thereby predicting other measurements. If a prediction obtained in this way 
fails, one would conclude that higher order effects have to be unusually large or that there 
is another deeper reason why a perturbative treatment does not apply. 

Such independent determinations of the same quantity of course probe the overall 
theoretical control that we have established. By themselves they do not tell us whether 
a failure found is due to unusally large higher order contributions or to a breakdown in 
duality. 

The fact that both the inclusive and exclusive methods for extracting |y(c6)| yield 
consistent values - and that the theoretical corrections one had to apply are both nontrivial 
and essential for the agreement - is such a test. We want to point to two other such tests 
that have become available, namely concerning the b quark mass and its kinetic energy 
expectation value. 

4.1 b quark mass 

The b quark mass has been extracted from beauty production at threshold in e~^e~ anni- 
hilation by several authors 0. Their findings can be stated in terms of two definitions of 
quark masses: 
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(i) The 'kinetic mass' is defined by 

1 = -77 Q + C' , (10) 

normahzed at 1 GeV; it is well-defined in full QCD and does not suffer from a 
renormalon ambiguity; equivalently one can use 



A(^) = hm^^^oo [M{Hq) - m§"(/.) 
where M{Hq) is the mass of the 0~ ground state. 



'Ill 



(ii) The pole or HQET mass which is a very popular choice, although it is not well-defined 
in full QCD since it suffers from the renormalon ambiguity. If appropriate care and 
caution are applied one can still use it in calculation; as a rule of thumb one has for 
its relatioship to the kinetic mass: 

Ahqet = A(l GeV) - 0.255 GeV (12) 

where the parameter A is defined in the same way as in (pT]): 

A^«^^ ^ hm^^^oo [M(i^Q) - m^°'^] (13) 

The results are completely consistent within the stated uncertainties of 1-2 % and can be 
summarized as follows: 

m^^"(l GeV)|,+e-^bb = 4.57 ± 0.06 GeV ^ A(l GeV) = 0.71 ± 0.06 GeV (14) 

or 

^pole ^ ^ g2 ^ Q Qg Qgy ^ ^HQET ^ g 45 ^ q gg Q^y ^^^-^ 

The techniques employed in the analysis differ somewhat from author to author; the full 
agreement in their findings is thus quite re-assuring. One should keep in mind, though, that 
these determinations share their experimental input to a large degree. The value stated 
in Eq.([T^) could thus be subject to some systematic shift from the true value. Arguments 
based on the small- velocity sum rules indeed suggest that m^*"'(l GeV) could lie a bit 
above 4.6 GeV. 

The b quark mass also affects the shape of lepton energy and hadronic mass spectra 
in semileptonic (and photon spectra in radiative) B decays. Its value can therefore be 
obtained from the measured lepton energy and hadronic mass moments, which encode the 
shape of these spectra. The DELPHI and CLEO collaborations have presented data as 
shown in the Figures ^ and Again it is pleasing to see that the different moments indeed 
yield completely consistent values although this is not truly surprising since they are highly 
correlated. DELPHI finds 

m^*"(l GeV)Uom = 4.59 ± 0.08 ± 0.01 GeV (16) 

^ A(l GeV) = 0.69 ±0.08 ±0.01 GeV 

mr'imom = 4.88 ±0.10 ±0.02 GeV (17) 
^ ^HQET ^ Q^4Q ^ Q^^Q ^ Q^Q2 GeV 
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M,(Mx) M,(Mx) M3(Mx) M,(E,) M,(E,) M3(E,) 




A (GeV) 



Figure 1: Contours in the A-Ai plane, from the DELPHI measurement of the first and 



second moments |TT|. The band defined by the almost flat curves in the upper part are 
from the second moment of the hadronic invariant mass, the steepest straight line is from 
the first moment of the hadronic invariant mass. The remaining curves are from the 
moments of the lepton energy, the steeper being from the first moment. The curve in the 
lower right is from the third moment of the hadronic in variant mass. 

It is again reassuring that their fit results are consistent with Eq.(|l^ CLEO measures 
truncated lepton energy moments and states their findings in terms of HQET parameters 

mr'lmom = 4.88 ± 0.03 ± 0.06 ± 0.12 GeV (18) 

^ J^HQET ^ Q gg ^ Q Q3 ^ Q gg ^ q -^^ GcV 

The main news here is how well the values for mb extracted from two sources agree that 
are quite different in their experimental as well as theoretical systematics, namely weak B 
decays and the electromagnetic production of beauty at threshold. 

4.2 Average kinetic energy 

A similarly pleasing picture emerges from determining the kinetic expectation value. Based 
on QCD sum rules and SV sum rules one had inferred for the infrared stable quantity 



10 





0.05 

0.10 

0.15 

0.20 

0.25 

^-0.30 

-0.35 

-0.40 

-0.45 
-0.50 



CM 
> 
CD 



1630802-010 



I IVW^I I ''J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

O*"^ Moment of 
Lepton Energy 

St 



1^* Moment 
_ of Lepton 
Energy 




la Theoretical 
Ellipse 



(Leptons) 
la Total 
Experimental 
Ellipse 



Moment of 
Photon Energy 
(b — sy) 



1 Hadronic _ 
Mass Moment 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
A (GeV) 



Figure 2: Contours in the A-Ai plane, from the CLEO measurements of the moments [12 . 



/i^(lGeV): 



filil GeV) ~ 0.45 ± 0.1 ( GeV)^ 



Using the rule of thumb for the relation to the HQET parameter Ai |T3 

- Ai ~ filil GeV) - 0.18( GeV)=^ 



(19) 



(20) 



this estimate translates into 

- Ai ~ 0.27±0.1(GeV)^ (21) 
The aforementioned DELPHI and CLEO analyses also yield values for this quantity, namely 



/i^(l GeV) = 0.31 ± 0.07 ± 0.02 ( GeV)^ DELPHI 
-Ai = 0.15 ± 0.07 ± 0.03 ( GeV)2 DELPHI 
-Ai = 0.25 ± 0.02 ± 0.05 ± 0.14 ( GeV)^ CLEO 



(22) 
(23) 
(24) 



The fact that the parameters extracted in different way and form different observables 
yield consistent values for the quark mass and the kinetic energy parameter indicates that 
no anomalously large higher order corrections or unexpectedly sizeable duality violating 
contributions are present. 
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5 Conclusions 



From all what we know currently from purely theretical considerations duality violations 
should be safely below one percent in the semileptonic branching ratio. This is likely to 
remain in the noise level of theoretical uncertainties due to terms of order 1 / ml and higher 
and of higher order perturbative contributions. Hence we do not see any need to assign 
some additional uncertainty to the extraction of Vcb from a possible duality violation in 
inclusive decays. This should not be seen as an ex cathedra statement. When more and 
more types of moments will be measured with more and more accuracy - even separately 
in the decays of B^, and B^ mesons additional constraints will be placed on the same 
set of heavy quark parameters. This will provide highly nontrivial tests of our theoretical 
control. 
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